Spontaneous symmetry breaking in cosmos: The hybrid symmetron 
as a dark energy switching device 
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We consider symmetron model in a generalized background with a hope to make it compatible 
with dark energy. We observe a "no go" theorem at least in case of a conformal coupling. Being 
convinced of symmetron incapability to be dark energy, we try to retain its role for spontaneous 
symmetry breaking and assign the role of dark energy either to standard quintessence or F(R) theory 
which are switched on by symmetron field in the symmetry broken phase. The scenario reduces 
to standard Einstein gravity in the high density region. After the phase transition generated by 
symmetron field, either the F(R) gravity or the standard quintessence are induced in the low density 
region, we demonstrate that local gravity constraints and other requirements are satisfied although 
the model could generate the late-time acceleration of Universe. 

PACS numbers: 98.80.Cq, 04.50.Kd, 04.50.-h 
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I. INTRODUCTION 

One of the inconsistencies of hot big bang model is 
associated with the observed primordial density fluctua- 
tion necessary for large scale structure in Universe. The 
paradigm which successfully addresses this problem along 
with the resolution of other theoretical issues such as flat- 
ness and horizon problems is inflation. The standard 
model of Universe is also plagued with age crisis which is 
late time inconsistency. In the standard model, the only 
known way to address this problem is provided by late 
time acceleration of Universe. This serves as an indirect 
evidence of the presence of an exotic matter of repulsive 
nature in Universe dubbed dark energy in the framework 
of standard lore. The discovery of the accelerating expan- 
sion of Universe 0, Q in 1998 supports this hypothesis. 
These observations tell us that the dark energy amounts 
to about 73% of the total energy budget of the Universe 
at present [3MTT| . 

The simplest cosmological model of dark energy is 
based upon cosmological constant (A) but the frame- 
work faces the most difficult conceptual and theoreti- 
cal issues associated with (A). With a hope to alle- 
viate some of these problems, a variety of scalar field 
models such as quintessence, phantoms, k-essence and 
tachyons have been investigated in the literature. p~2l — 

. Naively, the mass of slowly rolling quintessence field 
should be of the order Hubble rate in the present Uni- 
verse, Hq ~ 10 -33 eV. In case, it is much larger (smaller) 
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than Hq, it would mimic stiff matter (A). The incredible 
small mass scale of dark energy gives rise to formidable 
problems. In case, the scalar field is coupled to matter, as 
it should be in general, the propagation of the scalar field 
could generate the large correction to the Newton law and 
could be excluded by the local physics constraints. 

In order to address the problem related to local grav- 
ity constraints, three mechanisms of mass screening have 
been employed in cosmology. (1) Chameleon scenario: 
This mechanism operates with the field mass which be- 
comes dependent on the local density of environment 
such that the latter gets large in high density regime 
thereby leading to effectively decoupling of field from 
matter or suppression of fifth force [13, 13 • (2) Vain- 
shtein mass screening [19j : This mechanism is superior 
to chameleon and operates dynamically with non-linear 
derivative interactions. In the neighborhood of a massive 
body, the non-linear kinetic terms become strong leading 
to the decoupling of the field from the source in a large re- 
gion, around the massive body, specified by the so called 
Vainshtein radius. 

Recently, a very interesting, third screening mechanism 
closely related to chameleon dubbed symmetron is pro- 
posed in Refs. [2CH22l |. The idea of symmetron is related 
to the late time cosmic phase transition via spontaneous 
symmetry breaking. Similar to the chameleon model, 
the mass of the symmetron scalar field depends on the 
density of environment in a specific way due to its direct 
coupling to matter. For densities lower than some critical 
value, symmetron becomes tachyon and the symmetric 
vacuum state (<f> = 0) is no longer a true vacuum. In this 
case, there are two vacua with Zi symmetry which breaks 
soon after we choose one of these. The mass around the 
true minimum is well behaved. It is clear that if sym- 
metron is to be relevant to the dark energy, the phase 
transition should take place when the density of environ- 
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ment is low and the mass of the symmetron in the true 
vacuum is around H , which could be the first require- 
ment. The second requirement is that the symmetron 
should be invisible locally, that is, the fifth force induced 
by the symmetron should be negligibly small as compared 
to the Newtonian force of gravity. The second require- 
ment can be easily satisfied in the original symmetron 
model. As for the first, however, local physics imposes 
severe constraints on the symmetron mass which turns 
out to be quite heavy to derive late time acceleration. 
Although the symmetron might play some role during 
structure formation, the original idea of introduction of 
symmetron seems to be defeated. 

No doubt that the symmetron presents a beautiful 
idea and we believe that beauty cannot go for waste. 
Though the original idea of the symmetron does not seem 
to work for dark energy without unnatural fine tunings 
but the role of symmetron in cosmic symmetry break- 
ing could be retained. We propose a model where we 
use the symmetron to felicitate the cosmic phase tran- 
sition and the role of dark energy is played either by 
a second quintessence field or F(R) gravity which are 
switched on by the symmetron after the phase transition 
is over. In our proposal, the action reduces to that of 
the F(R) gravity dS-IH] (for reviews, see fllLIMI^) or 
standard quintessence field ip in the bulk, where the en- 
ergy density is very small but action reduces to the usual 
Einstcin-Hilbert one in the high density region, like in/on 
the earth, in the solar system, and in the galaxies, which 
makes the model consistent with the observations and 
local experiments. 

The plan of the paper is as follows. We, first, very 
briefly revisit the standard symmetron scenario, specially 
focusing on its failure for dark energy. As an attempt to 
to avoid the problem, we investigate a model where sym- 
metron couples to curvature scalar but observe a no go 
theorem for late time acceleration. Hence we propose a 
model where F{R) gravity couples with symmetron and 
show that there exists a parameter region which satis- 
fies possible constraints. The last section is devoted to 
summary and discussion. 



II. STANDARD SYMMETRON 



In this section, we briefly review the standard sym- 
metron model. This model is a generic attempt to imple- 
ment the original Ginsburg-Landau idea of phase transi- 
tions in cosmology. The model is based upon Higgs type 
scalar field potential with Z2 symmetry. In the high den- 
sity regime, the system resides in the symmetric ground 
state, symmetry breaks spontaneously in low density re- 
gion. The consistency of the model with local physics 
severely constraints the symmetron field to be relevant 
to dark energy [2l], [22j ■ 



A. Failure of symmetron to be dark energy 

Let us very briefly outline the basic features of the orig- 
inal symmetron scenario based upon the following Ein- 
stein frame UlUll 



S = 



V(<f>) 



(i) 



Here S m expresses the action of matter denoted by $ m . 
The scalar field equation due to its direct coupling to 
matter would involve density dependence 



- A3 dA 
P A '37, 



(2) 



where p is matter density in Jordan frame and is con- 
served in this frame. Since <j) couples conformaly to mat- 
ter, p is not conserved quantity in the Einstein frame. 
Usually a conserved quantity (p) in the Einstein frame 
is defined. The symmetron potential is invariant under 
Z2 symmetry (0 — > —<f>) and one preserves this symme- 
try in the effective potential necessary to implement the 
spontaneous symmetry breaking by making the following 
choice for A{4>) 
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2M 2 



where M is a mass scale in the model, 
potential then takes the following form 



(3) 



The effective 



(4) 



The mass of the field now depends upon the density of 
environment, naively, the field mass is given by, m 2 s = 
p/M 2 — p 2 . In high density regime, the mass depends 
upon density linearly, m 2 s ~ p/M 2 > 0. In this case, 
the system resides in the symmetric vacuum = 0. The 
requirement of local gravity constraints puts an upper 
bound on M and there is no a priori reason for it to 
be consistent with dark energy. In case of chameleon, 
there is more flexibility, the mass depends on density 
non-linearly. As shown in [HI, [HI, M < 10 _4 Mpi in case 
of symmetron. As the density redshifts with expansion 
and p falls below p 2 M 2 , tachyonic instability builds in 
the system and the symmetric state 0o = is no longer 
a true minimum; the true minima are then given by 



± 



p 2 - p/M* 



A 



(5) 



and the mass of the symmetron about the true minimum 
is, m s = \p2~p (at low density). Universe goes through 
a crucial transition when late time acceleration sets in 
around the redshift z ~ 1. One thus assumes that the 
phase transition or symmetry breaking takes place when 
p is around p cr , 
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HjM 2 , 
M 2 



m s 



M 



(6) 
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This means that m s > 10 4 Hq which is larger than the re- 
quired quintessence mass by four orders of magnitude. In 
this case, the field rolls too fast around the present epoch 
making itself untenable for cosmic acceleration. Invok- 
ing the more complicated potential with minimum with 
the required potential hight does not solve the problem; 
field goes oscillating for a long time and docs not settle 
in the minimum unless one arranges symmetry breaking 
very near to z = by unnatural fine tuning of param- 
eters [2l|, [22j ■ This also undermines the beauty of the 
underlying theory which is renormalizable and has an 
edge over chameleon theories which use potentials with 
complicated functional forms. In what follows we shall 
present a model which uses the A(/> 4 theory in a modified 
background. 



B. An effort to make symmetron compatible with 
dark energy that does not succeed 

Since the complicated choice of original potential will 
not help. We thus modify the effective potential by virtue 
of the modification of the background. We consider the 
following model, 
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F(cj>)R--{Vct>y-V{<t>) 
+S m [A 2 {4>)g^ 



(7) 



where F(<f>) is an arbitrary function of the field which 
should satisfy the obvious requirements of Z2 symme- 
try and the requirement that the model should not have 
instabilities. It is interesting to notice that if F = 
A 2 , the model ([7]) is equivalent to a minimally coupled 
quintessence field. 

Variation of action with respect to g^ v and <fi gives the 
following equations of motion 
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+ V^F - g^DF + -T^ , 
□0 =V(<t>)-RF'-^T. 



(8) 
(9) 



Taking the trace of Eq. ([8]) , we can express R through T 
and combinations of 4> and substituting the same in the 
field equation gives 



[F + 3F' 2 ) n<j> + F' 



3F" 



(V0) 2 

FV'-2F'V. (10) 



The second term in Eq. (|10[) is an extra kinetic term ob- 
tained due to field coupling with curvature. We should 
ensure that it docs not give rise to instabilities. For sim- 
plicity, it is better to switch it off but keeping modifica- 
tion of gravity alive, 



F' 3F" 



(11) 



In (fTTjl . an option F' = corresponds to the general 
relativity by Einstein, which is not desirable for us. The 
second option is interesting, 



ZF" + - = 



F(<f>) 



— ■ (12) 
12 v ; 



Since we like to have Zi symmetry to hold, we choose 
f3 = and for a, the obvious choice is, a = Mpj/2. We 
then have the expression for F(<p), 



F(4>) 



M 2 
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(13) 



We notice that F(<f>) > in generic range of <fi and cor- 
responds to the conformal coupling. The field equation 
with this choice becomes 



pi A' 

—F ) T + FV' — 2F'V . (14) 

2 A 



We now make use of generic functional forms of V((p) and 

A(4>), 



1 .2 i.2 , 1 \ /4 



V(^) = --^^ + -\^, A(<t>) = l + e-^ s . (15) 

The field equation (fT4"]) then takes the following form 
0<f> = dV c s/d(f> with 



V ee = 



1 

f 4 



P ~ - 

2 \M 2 

M 2 



&M 2 X 12M 2 M|j 



(16) 



In order to carry out the consistency check, it is instruc- 
tive to look at the expression of V e ff in low and high 
density regimes, 
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High density regime : V e g 
Low density regime : V c g 



e 

.2 ±2 



6M|] 



24M 2 M P ! 1 



(17) 
(18) 



Stability then demands that the effective potential should 
be bounded from below in any regime or the coefficient 
of </> 4 should be positive. Secondly the coefficient of 
the kinetic energy term should also be positive neces- 
sary to avoid the ghost instabilities. And these require- 
ments translate into restrictions on parameters, namely, 
e < (without no loss of generality we can take e = —1), 

1 -64>°- M>V6M Pl . 

Let us now consider spontaneous symmetry breaking 
which takes place when the coefficient of (f> 2 in Eq. (Tl7j]) 
turns negative which happens when 



Psb = P 



6M|] 



1 

M 2 



(19) 



Since we like to have the phase transition to take place at 
late times when Universe enters the regime of accelerated 
expansion, we shall identify psb with p CI 



Psb 



( 1 



\6M 



pi 



1 

M 2 



H 2 M 2 V 



(20) 

Since the effective mass of the field in the low density 
regime is m s = and we like to have it to be of the 
order of Hq, we define [i = aH where a = 0(1). Using 
the relation, M > \f&Mp\, we notice that 



a 2 H 2 



< a < 1/V6. 



(21) 

It is important to check whether or not the model is con- 
sistent with local physics or the local gravity constraint 
on symmetron mass are compatible with dark energy. In 
order to do that, we need to transform the action (0 into 
the Einstein frame, we have 
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jn E __ ( vvo 2 -y(vO 



+S m [A 2 (i;)g^,y m ] , (22) 



where the transformed potential and the coupling are 



given by 



A(i/>) = 



OM^Asinh 4 



o M pi , 
6e - tanh 



( H' 

/6 

i 

( 4' 



PI 



x cosh 



M 2 

f 4> 



pi 



(23) 



(24) 



and ip = Mpi\/6 axctanh ( ^ M 

For generic range of field values V <C Mp\, it can be 
approximated by 



^ 2 
2 V 



X 



4 y 2 



9M 2 y 



4 



^ 2 + J V 



Aty) ~1 + 



12Af2j 



M 2 
M 2 



=1 



12A/2j 



(25) 



(26) 



As we can see A is now a concave up function which con- 
tradicts the previously result e = — 1. Also if we push for- 
ward the analysis of the model, local gravity constraints 
impose a > 10 4 which do not match the cosmological 
condition derived previously (|2T|) . 

Let us try to understand why it happened. Broadly, 
the modified background brings two changes from the 
point of view of Einstein frame. First, the effective poten- 
tial, secondly, the coupling A 2 in the matter Lagrangian. 
Both terms are drastically modified but in the leading 
approximation, 4>/M <C 1 (4>/Mp\ ^ 1), we recover the 
standard symmetron action. Thus roughly, the mass 

scale M got replaced by (-1/M 2 + l/GM^)" 1 which 
also defines the mass of symmetron in the broken phase, 
see Eq. thereby leading to the same local gravity 

constraint as in the original symmetron model. Hence we 
see that the only difference would be a renormalization of 
the coefficients in the model, but the constraints will re- 
main same, and therefore the symmetron cannot play the 
role of a dark energy fluid. It is important to notice that 
we assumed tp <§C Mpi in order to derive (|26l25p . In fact 
the modified model proposed is equivalent to the origi- 
nal symmetron model because the field is slowly rolling, 
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hence A ~ 1. It is important to say few words about this 
generic feature of modified gravity models. 

Recently, a no go theorem related to the scope of mod- 
ified theories based upon chameleon/symmetron was dis- 
cussed in Ref. [32j |, We have not tried to avoid the "no 
go" theorem, our proposal is rather in agreement with 
Ref. [HI that late time acceleration is not the result of 
gravity modification but is driven by some quintessence 
field or a cosmological constant. In the subsection to fol- 
low, we add some clarifications on the important findings 
of Ref. [13 in the background of our proposal. 



Let us briefly check how it happens. We have the 
following translation between the Einstein and Jordan 
frames. 



a J (t J ) = A{4>)a E (t E ) , dt J = A^dt 1 



(27) 



and conformal time which is same in both the frames, 
dt = a(t)dr). In (f2"7|) . a J (a E ) is a scale factor and t J 
(t E ) is the cosmic time in the Jordan (Einstein) frame, 
respectively. 

Let us take the derivatives with respect to the Jordan 
cosmic time t J of a J (tj) = A(4>)a E (t E ) left right, 



III. SCOPE OF CHAMELEON/SYMMETRON 
SUPPORTED MODIFICATION OF GRAVITY 
FOR SELF ACCELERATION 

There are two ways of gravity modification: (1) Apart 
from the massless spin- 2 object, there is essentially at 
least an extra scalar degree of freedom which is exchanged 
apart from the graviton. In order to be relevant at large 
scales where modification is sought, the latter should 
gives rise to the effect of same strength as that gravita- 
tional interaction. This would then cause a havoc locally. 
The screening mechanism mentioned above address this 
problem by locally suppressing the exchange effects of 
the scalar degree of freedom. (2) The conformal coupling 
A{(f)) also modifies the strength of gravitational interac- 
tion. To pass the local tests, A(4>) should be very closely 
equal to one in high density regime in theories based 
upon chameleon/symmetron mechanism. As mentioned 
above, Universe has undergone a phase transition dur- 
ing < z < 1 which is a large scale phenomenon and 
one might think that the screening which is a local phe- 
nomenon should not impose severe constraints on how 
A(<j>) changes during the period acceleration sets in. It, 
however, turns out that the change the conformal cou- 
pling suffers as redshift changes from one to zero is negli- 
gibly small. Then a question arises, can such a conformal 
coupling be relevant to late time acceleration? It is well 
known that the de Sitter Universe is conformally equiv- 
alent to the Minkowski space-time. Has the conformal 
transformation changed physics? By 'physics', we mean 
the relationship between physical observables which is 
same in both the frames. In the Einstein frame we have 
the Minkowski space-time in which there is a scalar field 
sourced by the conformal coupling which is dynamical 
and couples to matter directly. The masses of all mate- 
rial particles are time dependent by virtue of A((j>). As 
a result, one would see the same relations between phys- 
ical observables in both the frames. The acceleration 
dubbed self acceleration is generic which can be removed 
(caused) by virtue of conformal coupling [321 ]. Simi- 
lar thing happ ens with Dvali-Gabadadze-Porrati (DGP) 
model [331 134| . Late time cosmic acceleration which can- 
not be affected by conformal coupling is caused by the 
presence of slowly rolling (coupled) quintessence and is 
not a generic effect of modified theory of gravity. 



d J (t J ) 




A TP A . TP 

T* a + A a 



(28) 



By time derivative "dot" ( ' ) of quantity in the Jordan 
(Einstein) frame, we mean time derivative with respect 
to the Jordan (Einstein) time t J (t E ). Multiplying this 
equation left right by a J = Aa E , we have, 



d J a J - d E a E 



-.E\ 2 , A E E 



The right hand side can be put in compact form by chang- 
ing the Einstein frame time to the conformal time from 
(d/dt E -> (l/a)(d/dr))) 




A r 



A{a E Y 



(a E Y 



(30) 



The first and last term are canceled with each other and 
the dependence on scale factor disappears: 



d J a J 



-EE 

a a = 



(A" 


A' 2 \ 




IT 


A 2 j 


-fir) 



(31) 



On the right hand side, "prime" ( ' ) denotes the deriva- 
tive with conformal time in the Einstein frame. Though 
the conformal time is same in both frames but coupling 
is defined only in the Einstein frame. 

Let notice that acceleration in the Einstein frame can- 
not be caused by conformal coupling, 



1 



6Mpj 



((P4 + 3P 4 ,)+p aP A(<j>)) 



(32) 



It is clear that in case acceleration takes place in the 
Einstein frame, it can only be caused by slowly rolling 
quintessence {p$ + 3P^ < 0). This implies that accel- 
eration in the Jordan frame and no acceleration in the 
Einstein frame is generic effect of conformal coupling or 
gravity modification. In this case, while passing from 
the Jordan to the Einstein frame, the acceleration is com- 
pletely removed, its affects in the Einstein frame are con- 
tained in the coupling such that the relationship between 
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physical observables are same in both the frames. The 
definition of self acceleration (32|, a E a 
then implies 



E - E < (d J a J > 0), 



(33) 



Since, A' = a J AA (AA - change over one Hubble (Jor- 
dan) time), it follows that 



,AA\ 
A ) 



3 « > a j ti J . 



AA 



1 dA\ 

ww)- (34) 



Integrating left right the above relation, we get [3/ 
AA 



A 



> 1 



(35) 



This quantity gives formally change of A over one Hubble 
time in the Jordan frame. As mentioned earlier, A is 
defined in the Einstein frame and naturally, the estimates 
of A from screening are obtained in that frame 



AA 
A 



/A 



1 



[ h e dt E ) I A 



(36) 



In case, AA is small or is of the order of one in the 
Einstein frame, it will be so in the Jordan frame. 

As demonstrated in Ref. (32[, screening imposes a se- 
vere constraint on the change of coupling during the last 
Hubble time, A A <C 1. Thus self acceleration cannot 
take place in this case. In most of the models supported 
by chameleon/symmetron screening, acceleration takes 
place in both frames such that d J a J and d E a E cancel 
each other with good accuracy or A A <C 1. In this 
case acceleration can only be caused by slowly rolling 
quintessence. In F(R) theories, the scalar field and the 
coupling both are made of F'(R) which imbibes the grav- 
ity modification. Since screening does not allow the con- 
formal coupling to felicitate self acceleration, the prob- 
lem simply reduces to coupled quintessence which one 
could deal with without really invoking F(R) theories. 
It is in this sense, the chameleon/symmetron supported 
modified gravity models have limited scope for late time 
cosmic acceleration. 



a complicated non-linear dependence on density of envi- 
ronment depending on the specific chameleon potential 
which allows chameleon to satisfy the local physics con- 
straints and also play the roll of dark energy. As for 
the symmetron modification, because the field is slowly 
varying in time, we have an approximate equivalence be- 
tween a modified symmetron action and the standard 
symmetron. Thus within the framework of (f> A theory 
with conformal coupling, the bound imposed by local 
physics on symmetron mass cannot be improved. How- 
ever, symmetron could still play the role of cosmic phase 
transition facilitator, the role of late time acceleration 
could be assigned to another field. 

In what follows, we shall consider a model based upon 
the following action, 

s=j dW^{^-i(v0) 2 -y(</>) 

dark energy 

}+S m [A 2 (4>)g^,y m ] . (37) 

where h{<f>) is a coupling function which satisfies the con- 
ditions 



h(0) = 0, h(<f> ) = l 



(38) 



for example, h{<j>) = (</>/^>o) m - 

We may rewrite the original symmetron potential in 
the following convenient form: 



V(</>) = A i 



(39) 



Here 



/j, 2 /AX. Notice that we added an effective con- 



stant (A0q) compared to the potential previously studied. 
As we will see, this constant will be always sub-dominant 
in cosmology and hence do not contribute to the acceler- 
ation of Universe. 

Let us note that in the high density regime such as 
the solar system or the galaxy, the action (f3"T|) reduces 
to Einstcin-Hilbcrt one with a sub-dominant symmetron 
field. On the other hand, in vacuum, the dark energy 
contribution to the action dominates. 



A. F(R) gravity switched by symmetron 



IV. THE HYBRID SYMMETRON MODEL 

The preceding discussion clearly shows that modifica- 
tions of symmetron action (within the framework of con- 
formal coupling at least) would not help to satisfy local 
gravity constraints and cosmological bounds suitable to 
late time evolution. 

Late us note that the symmetron mass linearly depends 
upon the density of the environment, the only parameter 
which remains to be constrained by local gravity tests 
in this scenario is the mass scale appearing in the cou- 
pling function. On the contrary, the chameleon mass has 



In this section, we propose a model where we use the 
symmetron to generate the cosmic phase transition and 
the role of dark energy is played by F(R) gravity which 
is switched on by the symmetron after the phase tran- 
sition. This mechanism could alleviate some fundamen- 
tal problems like the radiative corrections to the mass 
of the scalaron (or more generically, the mass of the 
quintessence field) [H], [29| and the Frolov singularity [3(| ■ 

The F(R) gravity is a famous model to explain the ac- 
celerating expansion of the present Universe. A problem 
in the F(R) gravity is the existence of the extra scalar 
mode which may give an observable correction to the 
Newton law. Such a correction could be observed by the 
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solar system test, the experiments on the earth, etc. In 
the present model, however, since the action reduces to 
that of the Einstein gravity in the region in solar sys- 
tems, in galaxy, or on the earth, the scalar mode does 
not appear and could be consistent with any local test of 
the general gravity. That is, the symmetron screens the 
extra scalar mode. 

On the other hand, in the vacuum, there appears the 
F(R) gravity, which generates accelerating expansion of 
Universe. A typical example is 

f(R) = A 2 R 2 . (40) 

Here A 2 is a constant. 

The ACDM model can be also realized by the F(R) 
gravity without introducing the real cold dark matter. In 
the model with a cosmological constant and the matter 
with the EoS parameter w, the scale factor behaves as 
a = a a e 9 ^ with 

*) = 5(T^) 1 "( Csi " 1 '(^ (i -'- ) ))' <41) 

Here, t s is a constant of the integration and C is a con- 
stant. The case w = corresponds to the ACDM model. 
The behavior of the scale factor a(t) can be realized by 
the F(R) gravity model by using the Gauss hypergeo- 
metric function |3l|. 

We should also note that the de Sitter space-time is an 
exact solution of the wide class of the F(R) gravity. 

In what follows, we investigate the profile of the scalar 
field (f> when the transition from the symmetric phase to 
the broken phase is induced; we explain how the local 
gravity constraints could be satisfied in the model under 
consideration. 

We may consider the symmetron like potential (|39l) . 
Then the field equation is given by 

□0 = 4A (0 2 - 0^ + ap0 - ti{<t>)F{R) , (42) 

where a = M~ 2 is the coupling to matter. We recover 
the previously studied model when F{R) = R and m = 2. 
In the following part, we do not consider this case. 

We know that F(R) ~ HfiM^ which impl ics a negli- 
gible contribution in the Klein Gordon equation. In fact 
when ~ 0, the F(R) term is trivially negligible, in 
the case where ~ 0o we need to impose the condition 
F(R) < acj>lp. We will define the range of viability for 
the parameters in order to satisfy this condition. Assum- 
ing that the conditions are satisfied, the equation reduces 
locally to the standard form without the F(i?)-term. 

We consider the static spherically symmetric solution 
in the flat space-time, the equation has the following 
form: 

hiift)-^-®** **' (43) 



We may assume 

(po when0^r<r o {u 
10 when r > Tq 

We also assume — 4A0q + ap > 0. 

First, we consider the behavior of in the region r <C 
tq and could be small. Then we may linearize (|43|) as 




The solution is given by 

Asinh (r^/— 4A0 2 , + ap ) 

4> = 0in = " L ■ (46) 

r 

Here A is a constant and we have assumed </> is finite at 
r = 0. 

Second, we consider the behavior of (f> in the region 
r 3> tq and </> ~ 4>q. By writing 

4> = 4>o + 64> , (47) 

we linearize the field equation (|4*5)) with respect to 5<j> as 
follows, 

Then the solution is given by 

4> = 4> out =0o . (49) 

r 

In the region r < ro, the solution PS]) grows up very 
rapidly when r increases. Inversely if we fixed a value of 
(f) for finite r < ro, 4> decreases very rapidly and goes to 
vanish as r decreases. 

Different from the original symmetry mod el dis cussed 
in Ref. [2l|, [22[ , we may assume the mass y/8\(f>Q of the 
symmetron in the bulk (low density regime) can be large. 
Then even in the region r > ro , (f> in the solution (|49[) goes 
to 0o very rapidly when r increases. The above behaviors 
tells us that <fi is almost constant except the small region 
r ~ ro- Then in order to estimate the constants A and 
B, we match the solutions </>;„ and O ut by imposing the 
following boundary conditions 

0inM = 0outM , 0LM = 0o ut M • (50) 

We find 
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.4 



B = 



yj 8A^g sinh (ry/— 4A0Q + c*Po) + V~ 4 ^o + a /°o cosh ^r-^/-4A(/)§ + ap 
(/) r e r \ / ^ (-^ sinh (?V-4A02 + apo ) + V /-4A<?!)2 + ap cosh (r v /-4A0g + apo 
Y/8A0Q sinh ^r v /-4A0Q + ap ) + \f-^4>l + ap cosh ^r v /-4A^Q + ap ) 



(51) 



r 



We may assume tq could be a radius of planet, star, or 
galaxy. Then ro should be much larger than the Comp- 
ton lengths 1 / y / 8A0f and 1 / yj — 4A0q 4- apQ. Then since 

sinh (ry/—AX(t)Q + apo^ ~ cosh (r\J— 4A0^ 4- apo) ~ 

e rV-4A^+ap / 2j we find 

20 o r o c-W- 4A ^+<Vo 



.4 



8A^ 



apa 



1 + 



(52) 



(53) 



We now consider the strength of the force generated 
by the scalar field <j). Since the strength of the coupling 
is given by a<j), we find 



Planck mass. Another constraint may come from the 
condition that the Compton length of the scalar field even 
in the vacuum should be much smaller than the size of 
galaxies re: 



V8A0f 



(57) 



We also require that h(4>)F(R) term does not affect the 
phase transition, 



(58) 



In case of the galaxy p ~ 10 5 H 2 Mp X and tq ~ 
(10- 7 - 10- 5 ) H Q \ 

Since H Q ~ 10~ 61 M P i and G = 8nM^, by introducing 
the new variables 



F^(r) = macj) 



dr 



(54) 



Here m is the mass of the particle receiving the force. In 
the region r > rg, by using (|49[) , we find 




On the other hand, the Newtonian force of gravity F g is 
given by F g {r) = GMm/r 2 , where M = 47rpo?~o/3. We 
notice that in the region specified by r 3> ro , the force F^ 
generated by the scalar field <fi can be neglected compared 
to the Newtonian force. 



B. Constraints on model parameters 

We now consider the constraints for the parameters. A 
constraint comes from the condition that the symmetry 
is restored when r < ro but broken in the vacuum: 



apo > 2A0Q > otp 



ccHqM^ . 



(56) 



Here p^ is the energy density of the vacuum, Hq is the 
Hubble constant in the present Universe, and Mpi is the 



x = 4>l , V = A0o , z 



1 

a 



yz < 10" 117 M^ , yz > VOT^M^ , 
y > lO" 10 ^/^ , xy > 10~ 122 M^ . 

By using, z = 10™Mp 1 , we obtain 



y > lCT 108 ^ , xy > 10 



PI ! 



-122 M 4 



PI 



We find xy > 10" 122 Af 4 , and for -14 < n < -9 
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-117-n i j-2 



Mp! > y > 10- luo M^i , x > 10 



(59) 



and putting ro = tq, the above constraints (|56j) . (J5TJ) , 
and (1581) can be rewritten as 



(60) 
(61) 



(62) 
(63) 



-5+n M 2 



PI I 

(64) 



and for n < —14 



10-iir-nj^ > y > lo-iM-njiflj , a: > 10~ 5+n Af ^ . 



(65) 



We therefore confirm that we have always a range of vi- 
ability of the parameters, where all the constraints are 
satisfied. 
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C. Phase transition in the early Universe 



and for n < — 14 



As mentioned earlier, the phase transition should have 
occurred in rather early Universe, not in the late Uni- 
verse. We now investigate how the transition could have 
occurred and if there could be any problem or not. The 
potential ([3"9")) tells that the critical density p cr , where the 
phase transition occurs, is given by 



4A</> 2 



4yz. 



(66) 



When the phase transition could have occurred, the Hub- 
ble rate H CT is given by 



H 2 - K n yZ 
H a - ~^P" 



TUT 2 



(67) 



Here we have used the notation defined in (|59[) . As we 
are interested in the behavior when <j) ~ <f>Q, by using 
(|47p. we linearize the scalar field equation as 



(68) 



dt 2 ■■3^ + 8A^ = 0. 



Equation (|68[) has a form analogous to the equation of 
motion of the harmonic oscillator with drag (air resis- 
tance) : 



fx 
It 2 



27— + w 2 x = 0. 
at 



(69) 



Here 7 and uj are positive constants and x expresses the 
position of the harmonic oscillator. As well-known, when 
7 2 > uj 2 , the amplitude |x| decreases without oscillation 
and when j 2 < uj 2 , the amplitude \x\ decreases with os- 
cillation by the angular velocity y 1 'uj 2 — ■y 2 . Hence the 
time scale relevant to the decrease of amplitude is given 
by T doc = I/7. Then Eq. ([68]) tells us that S(j) would 
vanish without oscillation if 9H 2 = > 8A0 2 , = 8y, 

that is, -rrr > 1. On the other hand, Sip vanishes with 

oscillation if j^- < 1. By writing z = 10"A/pj as before, 

and since n < —9, we find 8<j) vanishes with oscillation. 
The time scale of decreasing the amplitude is given by 

T dcc = l/H~M Ply /yZ. (70) 

Then Eq. (|64l65p tells when -14 < n < -9 

lO 58 ^ 1 < T dcc < lO 54 ^"/ 2 ^ 1 , (71) 

and for n < — 14 

lCPil/pj 1 < T dcc < lO^il/p/ 1 . (72) 

Since Mpj 1 ~ 10 -44 sec ~ 10~ 51 years, we have for — 14 < 

n < -9 



10 7 years < T dcc < I0 3 ~ n / 2 years , 



10 7 years < T dcc < 10 10 years . 



(74) 



We therefore conclude that it could take 10 million as 
minimum or 10 billion years as maximum for decreasing 
the oscill ation s. The period Tp of the oscillation is given 
by 27T/V8A0 2 - 1/jy. Therefore Eqs. ([751731) gives 



10 7+ 2 years < T P < 10 3 years 
10 7+ * years < T P < 10 10+ * 



years . 



(75) 
(76) 



for — 14 < n < —9 and n < — 14 respectively. For 
instance, if we choose n = —14, we have 1 year < Tp < 
10 3 years. Therefore the oscillations are very slow and do 
not give rise to particle production. 

Although the phase transition did not occur during 
the accelerating expansion, that is, when the rcdshift is 
z ~ 1, we may estimate when it could have occurred. 
Since we assume po ~ 10 5 H^Mp, > p CI > and 



(73) 



p cx a , we may assume the phase transition could have 
occurred when the redshift z — — l 10 (ao is the 
present value of the scale factor a), which corresponds 
to a few million years after the Big-Bang; the redshift of 
the cosmic microwave background radiation corresponds 
to z - 1000. 

In the early universe, the Zi symmetry could be re- 
stored. Then Eq. ([3"9"[) tells us that V{<j>) gives a contri- 
bution to the energy in the vacuum by A0 4 .. The con- 
tribution is, however, not so large. Equation ([59"[) tells 
A</>q = xy. Using the previous results, we find A</>q = 
xy > 10~ 122 (Mpi) 4 . On the other hand, for an example, 
the energy density of the matter in the present universe 
is 10~ 124 (Mpj) 4 . When z > 10, since the matter den- 
sity could be 10- 124 x (1 + zf (M P1 ) 4 - 10~ 121 (A/ P1 ) 4 
Then the magnitude of V(0) = A0 4 , is comparable with 
the energy density of matter but could not be dominant. 
Therefore the contribution from the constant A$>q = xy 
can always be neglected. 

Let us note that £ dar k energy in ([37)) could also be rep- 
resented by a regular quintessence Lagrangian C v . The 
dynamics during phase transition then proceed in the fol- 
lowing manner. In the beginning as <f> begins to evolve 
from the symmetric state <f> = towards the symmetry 
breaking vacuum, the coupling function h(<j>) is negligi- 
bly small and as a result, the energy density p v ~ 0; the 
major contribution to the Hubble parameter is made by 
p,!,. Hence the field ip is frozen due to the Hubble damp- 
ing. As 4> evolves towards the minimum and h(<f>) builds 
up making the field ip alive which then begins rolling. 
Soon after <f> reaches a critical value (j> CI , the energy den- 
sity in ip begins to dominate and soon overwhelms the 
symmetron field which then freezes and its role ceases in 
the cosmic dynamics thereafter. The field ip then derives 
late time acceleration of Universe. 
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V. SUMMARY AND DISCUSSION 

In this paper, we examined the modified symmetron 
models with an aim to reconcile the latter with late time 
cosmic acceleration. We have investigated a symmetron 
type framework coupled to curvature scalar and to mat- 
ter. In this case the effective symmetron potential drasti- 
cally differs from the one in original scenario. We exam- 
ined the consistency of the model with local physics and 
found constraints on the symmetron mass similar to the 
one in the standard model. The latter is related to the 
fact that local gravity tests are insensitive to the details 
of the effective symmetron potential and the new mass 
scale appearing in the over all coupling in the leading 
approximation (<fi/M -C 1 , 4>/Mp\ <C 1) in the model 
under consideration is same as the one that defines the 
symmetron mass. We therefore conclude that there is a 
"no go" theorem for dark energy symmetron as long as 
we confine to conformal coupling, situation might (might 
not) change in the disformal case. 

Being inspired by the beauty of cosmic symmetry 
breaking, we have made a proposal in which symmetron 
felicitates the phase transition in low matter density 
regime ; the symmetry is restored in high density ensur- 
ing the compliance of the model with local gravity con- 
straints. The role of symmetron field ceases after phase 
transition thanks to its coupling to an F(R) gravity ac- 
tion or quintessence field. 

In case of F(R) symmetron or cquivalently a 
quintessence field, the model reduces to the standard 
Einstein gravity in the high density region whereas F(R) 
gravity (or quintessence) is switched on in the low den- 
sity region; the transition is generated by the symmetron 
field. Since in our proposal, symmetron field is responsi- 
ble for phase transition only, its mass could be large even 
in the vacuum and therefore the local gravity constraint 



can easily be satisfied, which is different from the origi- 
nal symmetron model pll . l22j . We have also investigated 
other constraints and we have shown that there exists a 
parameter region which satisfies all constraints. In the 
model under consideration, the transition could have oc- 
curred around z ~ 10 - the dark age. After the transi- 
tion, symmetron field oscillates around the true minimum 
but the time scale of decreasing of amplitude varies be- 
tween 10 million to 10 billion years which is much smaller 
compared to the Hubble scale. On the other hand, the 
period of the oscillations varies from one to one thousand 
years, which could be large enough from thermodynami- 
cal point of view. Thus the oscillation would not conflict 
with the observed cosmology. 

Last but not least, it might be interesting to enlarge 
the symmetron framework to disformal type of set up 
to check whether the disformal symmetron can reconcile 
with dark energy. 
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